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ABSTRACT 

We consider a 'brane-world scenario' recently introduced by Dvali, Gabadadze 
and Porrati, and subsequently proposed as an alternative to a cosmological con- 
stant in explaining the current acceleration of the universe. We show that, con- 
trary to these claims, this particular proposal is already strongly disfavoured by 
the available Type la Supernovae, Cosmic Microwave Background and cluster 
data. 

Subject headings: Cosmological parameters — cosmology: theory — methods: data 
analysis — supernovae: general 

1. Introduction 

At a time when observational cosmologists are finally pinning down some crucial cosmo- 
logical parameters (Perlmutter et al. 1997, 1999; Riess et al. 1998; de Bernardis et al. 2000; 
Lange et al. 2001; Netterfield et al. 2001; Hanany et al. 2000; Balbi et al. 2000; Stompor et al. 
2001; Halverson et al. 2001; Pryke et al. 2001; Percival et al. 2001), theoretical cosmologists 
have increased and diversified their efforts to try to provide some more solid connections 
between particle physics and cosmology. 
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In this context, a topic of much recent interest has been that of the so-called 'brane- world 
scenarios' (Arkani-Hamed, Dimopoulos & Dvali 1998; Randall & Sundrum 1999; Binetruy et 
al. 2000; Maartens 2001; Rubakov 2001). The general principle behind such models is that 
the ordinary particles live on a three-dimensional surface (commonly called a 3-brane, or 
simply 'the brane'), which is embedded in a larger space ('the bulk', which may or may not 
be compact and might even have an infinite volume) on which gravity can propagate. An 
observer on the brane will measure four-dimensional gravity up to some corrections which, 
given the weakness of gravity, can in general be made small enough not to conflict with 
observations without tweaking with model parameters too much. 

At present the topic is young enough that the main drive is still to try to explore all 
remotely viable model-building possibilities without worrying too much about the conse- 
quences. However, some of the proposed models are already developed enough that they 
can start to be put through the sieve of specific cosmological observations. In this paper we 
will provide what we believe to be the first detailed analysis of this kind for a brane world 
scenario. 

We will consider a particular solution of a brane world scenario originally introduced by 
Dvali, Gabadadze & Porrati (2000), and further studied in (Deffayet 2001; Deffayet, Dvali & 
Gabadadze 2001) — we shall refer to it as the DGP model for simplicity. It's a five-dimensional 
brane- world model with a non-compact, infinite- volume extra dimension. The usual four- 
dimensional gravity is recovered on the brane for scales below a 'characteristic radius' r c , due 
to a four-dimensional Ricci scalar being induced on the brane. However, at larger scales this 
becomes sub-dominant, and one will effectively see five- dimensional Einstein-Hilbert gravity. 

The above effect can obviously have dramatic cosmological implications. A particu- 
larly interesting solution was first found by Deffayet (2001) — and then generalized by Dick 
(2001) — and then further studied in (Deffayet, Dvali & Gabadadze 2001). It describes a uni- 
verse which at late times is accelerated on scales larger than r c . This is an effect of the bulk 
gravity, in the sense that observers on the brane will see no cosmological constant. Hence 
this is another interesting alternative way to explain the current acceleration of the universe, 
which is strongly indicated by Type la supernova observations (Perlmutter et al. 1997, 1999; 
Riess et al. 1998), without resorting to a cosmological constant — for earlier alternative ex- 
planations, see (Dev, Sethi & Lohiya 2000; Behnke et al. 2001; Mannheim 2001). Note that, 
unlike most other known brane worlds scenarios, here the early evolution of the universe is 
the standard one while the late evolution is different. Also, unlike other alternative theories 
of gravity (introduced in other contexts) here gravity will become weaker on large enough 
scales. These two points will be important in what follows. 

In (Deffayet, Dvali & Gabadadze 2001) the authors argue that this is an intrinsically 
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higher-dimensional effect, at least in the sense that one can not mimic it with arbitrary high- 
derivative terms in ordinary four-dimensional gravity. This turns out to be both a blessing 
and a curse, for on one hand it means that one can extract quite distinctive observational 
predictions, but on the other hand it also implies that it's quite easy to rule it out. In (Def- 
fayet, Dvali & Gabadadze 2001) the authors claim that the model's alternative explanation 
for the current acceleration of the universe agrees with all existing cosmological observations 
(or, more accurately, that it is currently indistinguishable from the standard scenario). In 
what follows we shall show that this is not the case. Indeed, currently existing data is al- 
ready sufficient to make this alternative explanation for acceleration strongly disfavoured 
when compared to the standard one. 

The plan of the paper is as follows. In Sect. 2 we provide a brief summary of the features 
of the DGP model which are relevant for our discussion. We then proceed to analyse the 
accelerating solution in the light of the Type la supernovae data in Sect. 3. In Sect. 4 
we cross-check the results of this analysis with other cosmological datasets, and finally we 
present our conclusions in Sect. 5. 

2. The Model 

Here we briefly describe the brane-world scenario introduced by Dvali, Gabadadze & 
Porrati (2000), and further studied in (Deffayet 2001; Deffayet, Dvali & Gabadadze 2001)- 
we shall henceforth refer to it as the DGP model. Our discussion will be somewhat simplified, 
as we will only focus on the features that are relevant for our subsequent analysis — the reader 
is encouraged to consult the original references for a more detailed discussion. 

Our three-brane is embedded in a five-dimensional spacetime with a non-compact, 
infinite-volume extra dimension. Particles in the standard model are confined to the brane, 
and brane fluctuations are neglected. There is essentially one free parameter in the model, 
which is the 'five- dimensional Planck mass', denoted M 5 . Note that one must assume that 
the standard model cut-off doesn't coincide with M 5 — in fact, it must be much larger, so 
that the physical interpretation of M5 is not quite trivial. 

The four-dimensional Planck mass will be denoted M 4 , and is related to the usual grav- 
itational constant through 87rG 4 = M^ 2 . Unlike in other brane world scenarios, here the 
two masses M 4 and M 5 need not be related. We note that there is a somewhat technical 
problem with the model (Deffayet 2001) which implies that if one defines Newton's con- 
stant via a standard Cavendish-like experiment, then the so-defined Gn doesn't necessarily 
coincide with G 4 . This would obviously contradict standard tests of General Relativity, as 
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was already pointed out in (Deffayet 2001). Possible ways to circumvent this problem have 
been claimed (Dvali, Gabadadze & Porrati 2000; Deffayet, Dvali & Gabadadze 2001). In 
any case, we shall neglect this aspect in what follows, since our present purpose is to discuss 
'cosmological' (as opposed to 'local') tests of the model. 

The usual four-dimensional gravity is recovered on the brane for scales below a 'charac- 
teristic radius' r c , which is given by 
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This is due to a four-dimensional Ricci scalar being induced on the brane. However, at 
larger scales this term becomes sub-dominant, and one will effectively see five- dimensional 
Einstein-Hilbert gravity. Therefore gravity becomes weaker on large enough scales. This is 
to be contrasted with models where one modifies gravity on large scales in order to solve, for 
example, the dark matter problem: in that context, one requires stronger gravity on large 
scales. From this it immediately follows that one can impose a simple constraint on M5, 
since the characteristic radius must at least be as large as the present Hubble radius. 

Cosmological solutions in this model were first studied in (Deffayet 2001) — but see also 
(Dick 2001). One finds that the Friedmann equation on the brane has the following form 
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where k — 0, ±1 is the spatial curvature and e = ±1 corresponds to two different brane 
embeddings in the bulk spacetime. On the other hand, the energy conservation equation has 
the standard form, 

f t +3H(p + p) = 0. (3) 

At early enough times the density term dominates the Friedmann equation, and hence 
one obtains (at least to first order) the standard cosmological evolution, namely 3H 2 = 
RitGip. The additional bulk-induced term will become important when H~ l ~ r c . Then 
the subsequent evolution depends on the sign of the parameter e. In the e = — 1 branch the 
universe switches into a full five- dimensional gravity regime, where the Friedmann equation 
looks like H oc p — something that is typical of many brane world scenarios. On the other 
hand, in the e = +1 branch something rather more interesting happens. There is a 'self- 
inflationary' solution with H ~ r~ x . What happens is that the additional curvature term 
on the brane appears as a source for bulk gravity, and thus can cause acceleration on the 
brane. In other words, an observer on the brane will see the universe being accelerated on 
scales larger than r c . 



- 5 - 



Note, however, that this solution does not require any other energy source on the brane— 
so in this sense this is indeed a higher-dimensional effect. In particular, no cosmo logical 
constant is needed on the brane, so this is an interesting alternative way to explain the current 
acceleration of the universe, which is strongly indicated by Type la supernova observations 
(Perlmutter et al. 1997, 1999; Riess et al. 1998), without resorting to a cosmological constant. 
A simple, 'back-of-the-envelope' constraint comes from the fact that we want the universe 
to be at this crossover stage at about the present epoch if the alternative proposal for the 
acceleration of the universe is to be viable, hence Hq 1 ~ r c . This then naively implies that 
the five-dimensional Planck mass should be of the order of 

M 5 ~ 10 - lOOMeV . (4) 

We finally note that high-energy processes place almost no constraints on this mass 
scale M 5 , basically because up to about the present epoch the universe evolves as normal. 
Indeed the only constraint comes from the measurement of the Newtonian force, which only 
implies the very mild 

M 5 > l(T 3 e\/ . (5) 

As we shall see in the following section, much more stringent constraints can be derived 
using cosmological observations, if one assumes that the accelerating solution is valid. We 
will start by an analysis of the Type la supernovae data, which is described below. We will 
then contrast the results of this analysis with other cosmological constraints. 



3. Supernovae Data Analysis 

We begin by evaluating the luminosity distance as a function of the cosmological pa- 
rameters for our model. The Friedmann equation (2) can be rewritten as a function of the 
red-shift 1 + z = a /a to give: 



HI 




where 



respectively represent the fractional contribution of curvature, the bulk-induced term and 
the other components in the Friedmann equation. In equation (6) the sum is over all the 
components of the cosmic fluid with an equation of state p a = w a p a . From now on we 
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will consider only one component of non-relativistic particles together with the bulk-induced 
term, in which case equation (6) becomes 



HI - 



k (l + z) 2 + (v^ + v /fi r C + ^ m (i + ^) 3 ) 2 ; (8) 



in particular, at the present day one must have 

(9) 

For a flat universe Qk = and so in this case the two cosmological parameters are related 
by (9) 

1 _ o x 2 



It is also possible to show that closed universes with 

\Q k \ 3/2 >8Q m QH 2 (11) 

do not have a big bang. These universes avoid the big bang singularity by bouncing in the 
past. We shall disregard such universes in most of what follows. Another useful benchmark 
is the redshift at which the universe switches from deceleration to acceleration, or in other 
words the redshift for which the deceleration parameter vanishes. For a flat universe Qk — 
it's easy to show that the following exact result holds 

m \ 1/3 

note that in this (flat) case Vt rc and Q m are not independent parameters — they are related by 
(10). For comparison, the result in the standard case, again taking a flat model (f2 m +f2 A = 1) 
is 

< 1+ '>~=S:f ! (13) 

we shall return to these quantities in the following section. 

It is also straightforward to show that in a Friedmann- Robert son Walker (FRW) universe 
the luminosity distance is given by: 

d L = Ho\i + ^r i/2 Mi^r /2 ^) , (14) 

where 



dc 



=H °L'm- (15) 
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and Sk is sinh if Q k > and sin if Q k < 0. Note that the assumption of a FRW universe 
is the only one needed to derive (14). Specifically, we do not need to assume the validity 
of General Relativity or to specify a priori the detailed contents of the universe. This is of 
course crucial in our case. 

We estimated the cosmological parameters using the combined data of two independent 
teams — thus making up a dataset of 92 different supernovae — using the procedure described 
in (Wang et al. 2000; Wang & Garnavich 2001). The measured distance modulus for a SN 
la is 

$ = 1$ + *® (16) 

where $ is the theoretical prediction 

(17) 

and is the uncertainty in the measurement, including observational errors and intrinsic 
scatters in the SN la absolute magnitudes. 

We denote the parameters to be fitted as s and estimate them using a x 2 statistic, with 
(Riess et al. 1998) 



flp [Zi\S) - fl Q 



x 2 W = E i ? -rj ~ 5 £ j 3 ' < 18 > 



where cr w is the estimated measurement error of the distance modulus, and a mz is the 
dispersion in the distance modulus due to the dispersion in galaxy redshift, a z , due to 
peculiar velocities and uncertainty in the galaxy redshift. The probability density function 
(PDF) for the parameters s is 

p(s) oc exp f-y) • ( 19 ) 

The normalized PDF is obtained by dividing the above expression by its sum over all possible 
values of the parameters s. In the particular case of our model the cosmological parameters 
are H Q , Q m and Q rc . The probability distribution function for the parameters fl m and Q rc is 
obtained by integrating over all possible values of H , and the results are displayed in Fig. 
1. 

We have also studied the improvements in the parameter estimation using supernovae 
which are expected from future studies of cosmic acceleration. Following (Weller & Albrecht 
2001) we assumed a future dataset similar to one proposed for the SNAP satellite 4 . This has 



4 SNAP home page at http://snap.lbl.gov 
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the magnitudes of 50, 1800, 50 and 15 supernovae in the red-shift ranges from z = — 0.2, 
z = 0.2 — 1.2, z — 1.2 — 1.4 and z — 1.4 — 1.7 respectively. The statistical error in magnitude 
is assumed to be a = 0.15 including both the estimated measurement error of the distance 
modulus and the dispersion in the distance modulus due to the dispersion in galaxy redshift. 
The supernovae dataset was generated assuming that we live in a standard FRW universe 
with cosmological parameters Q m = 0.3 and Q\ = 0.7. Fig. 2 shows the corresponding 
results. 



4. Results and Discussion 

A number of interesting features are apparent from Fig. 1. Firstly, the likelihood 
analysis of the supernovae data is degenerate in the fl rc — Q m plane, approximately following 
a line of the form 

vt rc ~ -n m + — . (20) 

This is to be compared with the standard cosmological scenario, where the degeneracy is 
approximately along 

Ov ~ + £ ■ (21) 

Hence one can say that for any given value of Q m , the value of Q rc which provides the best- 
fit to the supernova data is always lower than the corresponding value of Reversing 
the argument, one could also say that for a given value of the density of the accelerating 
component (a cosmological constant in the standard case or the bulk-induced term in the 
DGP model) the DGP model requires a higher matter density in order to fit the supernova 
data. Note that the two and three sigma likelihood contours are quite close to each other, 
and relatively distant from the one sigma contour. This indicates that with the currently 
available data there is an elongated 'degenerate best-fit plateau', and beyond this plateau 
the likelihood drops quite abruptly. 

In any case, just by looking at Fig. 1 one might think that there is a rather comfortable 
range of matter densities which would give models in agreement with observation. However 
this is not the case as there are other cosmological constraints that must be met. In partic- 
ular, the most recent CMB data (de Bernardis et al. 2000; Lange et al. 2001; Netterfield et 
al. 2001; Hanany et al. 2000; Balbi et al. 2000; Stompor et al. 2001; Halverson et al. 2001; 
Pryke et al. 2001) gives a strong indication that the universe is spatially flat or very nearly 
so. The current constraint is 

Q tot = 1.00 ±0.05; (22) 
note that only fairly weak priors are needed to derive this constraint (refer to the original 
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CMB papers for the analysis details). Combining this with the supernova analysis this leaves 
a much smaller range of allowed models. At the 68% confidence level, the allowed range of 
matter densities is approximately 

tt m = 0.20 ±0.05, (23) 
while at the 99% confidence level it is 

0„ = 0.2 ±0.1. (24) 

Note that this result is quite robust — for example, we have checked that it is unchanged if 
the likelihood analysis is restricted ab initio to flat universes. 

The final piece of observational evidence that we shall use are dynamical measurements 
of the total mass density — see Turner (2000a,b) for a discussion of the state-of-the-art. In 
particular the ratio of baryons to the total mass in clusters has been determined using both 
X-ray measurements (Mohr et al. 1999) and Sunyaev-Zel'dovich measurements (Carlstrom 
et al. 1999). One respectively obtains 

f x = (0.075 ± 0.007)/T 3/2 (25) 

from X-ray measurements, and 

f 3Z = (0.079 ±0.010)/T 1 (26) 

from the Sunyaev-Zel'dovich effect. In both cases, various dozens of sources have been used 
in the analysis. If one assumes that clusters are a fair sample of the matter content of the 
universe (which is very reasonable given their large size) and uses the latest value of the 
baryon density at nucleosynthesis (Buries, Nollett & Turner 2000) 

tt b = (0.020 ±0.002)/T 2 , (27) 

together with the value of the (re-scaled) Hubble constant h obtained by the HST Key 
Project (Freedman et al. 2000), 

h = 0.72 ±0.08, (28) 
one finally obtains the (rather conservative) estimate 

n m = 0.35 ±0.07. (29) 

More recently, there have been claims of an even narrower (though perhaps slightly opti- 
mistic) range (Turner 2001) at the one sigma level 



Vl m = 0.330 ±0.035. 



(30) 
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Note also that there are various other sources of supporting evidence that are consistent 
with the above value, including studies of the evolution of cluster abundances with redshift, 
measurements of the power spectrum of large-scale structure (such as the recent preliminary 
2dF results (Percival et al. 2001)), analyses of measured peculiar velocities as they relate 
to the observed matter distribution, and observations of the outflow of material from voids. 
A discussion of the assumptions and techniques of each method can be found in Turner 
(2000a,b). 

Hence the DGP model's proposal for the acceleration of the universe requires a value of 
the matter density that is inconsistent, at least at the two sigma level, with the observation- 
ally estimated matter density of the universe. Together with the fact that gravity becomes 
weaker on large enough scales, this presents a serious problem. Note that if the required 
mass was larger than the standard case, one could perhaps argue that there was some matter 
in an yet undetected form. Indeed, the fact that in the DGP model gravity becomes weaker 
on large enough scales could then be used to obtain a relatively simple explanation. How- 
ever, since the observationally acceptable range of masses is lower than the standard model, 
no explanation of this kind is possible. In any case, one should recall that the evolution of 
the universe should be as standard in the DGP model up to very recent times, eg in what 
concerns the Friedmann equation for example — which also places strong constraints on any 
attempts to 'get rid of some of the matter. 

And finally, there is yet another hurdle for this model to overcome. In (12) we derived 
the redshift at which the universe switches from deceleration to acceleration, for the case of 
a flat universe. We plot this redshift, for the range of matter densities given by (24), in Fig. 
3, together with the analogous curve for the standard model. The problem is now apparent: 
for the specified range of matter densities, the redshift of turnaround decreases as the matter 
density increases. On the other hand, for a given value of the turnaround redshift, the 
required matter density is always lower with a bulk-induced term replacing a cosmological 
constant (the solid curve in Fig. 3) than in the standard model with a cosmological constant 
(the dashed curve). Now, the latest supernova data (Perlmutter et al. 1997, 1999; Riess et al. 
1998, 2001; Turner & Riess 2001; Avelino, de Carvalho & Martins 2001) indicates that the 
universe switched from deceleration to acceleration at a redshift in the interval (at one-sigma 
level) 

0.6 < ^ =0 < 1.7. (31) 

Note that for values of the matter density close to the upper limit Q m ~ 0.3 the predicted 
redshift of turnaround is already smaller than this range. 

This, therefore, is the dilemma of these models. On the assumption of a flat universe, 
a very low matter density is needed so that acceleration starts early enough. This is in fact 
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confirmed by the simulation of the supernova analysis for a SNAP-class dataset, which is 
shown in Fig. 2: closed models are still favoured, though flat ones are still possible at around 
the two sigma level. However the range of possible matter densities is significantly reduced. 
Note that in generating the SNAP dataset we have assumed a standard FRW universe with 
Q m = 0.3 and Q\ = 0.7. Hence, according to our discussion above, if we fit that dataset to 
the accelerating DGP model then the preferred value of the matter density will come out 
lower. This is a trivial consequence of the fact that the type la supernova analysis method 
is basically a cosmological 'accelerometer'. The point is that, even with the data available 
today, such low values are already strongly disfavoured by dynamical measurements of the 
total mass density in the universe. 

On the other hand, even if one would be willing to admit that such values were allowed 
on the grounds of dynamical measurements alone (implying a much smaller value of dark 
matter than in the standard model), they are expected to run into serious difficulties when 
it comes to density fluctuation growth and the evolution of large-scale structures (which can 
now be probed much beyond Mpc scales both by direct surveys and through gravitational 
lensing), again because of the weaker gravity on large enough scales. For the DGP scenario 
to be viable the characteristic scale would be of order r c ~ Hq 1 , but obviously the effects of 
weaker gravity would be felt on smaller scales than this. Indeed this point has already been 
made on rather general grounds (though only for the case of sub-horizon modes) by Uzan & 
Bernardeau (2000), and we shall return to it in more detail elsewhere. 

5. Conclusions 

In this paper we have considered the cosmological consequences of the brane world 
model of Dvali, Gabadadze & Porrati (2000), and its proposed alternative explanation for 
the current acceleration of the universe. We have shown that, contrary to recent claims 
(Deffayet, Dvali & Gabadadze 2001), this proposal is already strongly disfavoured by existing 
cosmological datasets, at least at the two sigma level. In order to be consistent with CMB 
and supernova data one would need a very low matter density Vt m ~ 0.2. Even if this was 
allowed by dynamical measurements (such as cluster data), such a low density ( and hence 
such a small amount of dark matter) together with the fact that gravity is weaker on large 
enough scales would make it difficult to produce a consistent structure formation scenario. 

The lesson to be learned from this exercise is twofold. Firstly, no matter how interesting 
or mathematically clever one's favourite particle physics model of the universe might be, the 
first hurdle towards credibility consists in deriving falsifiable cosmological predictions from it. 
And secondly, the currently available cosmological observations are already powerful enough 
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to impose tight constraints on a wide range of possible models, especially when various 
cosmological datasets are combined — which is a sign that the era of precision cosmology has 
indeed started. We hope that other brane world scenarios can be brought into the realm of 
cosmological testability in the near future. 
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Fig. 1. — The probability distribution function for the parameters Vt m and Vt rc in the DGP 
model, for the presently available dataset of 92 Type la supernovae — see the text for detailed 
description of the method. The 68%, 95% and 99% confidence contours in the Vl rc — Q m plane 
are shown, as well as the line separating closed universes with big bang from no-big-bang 
'bouncing' ones and the one separating closed and open universes (ie, denoting the flat ones). 
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Fig. 2. — A simulation of the analysis with a SNAP-class dataset generated assuming that we 
live in a standard FRW universe with cosmological parameters Q m = 0.3 and fl\ = 0.7 — see 
the text for the other assumptions involved. Contours and boundary lines are as in Fig. 1. 
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Fig. 3. — The redshift of 'turnaround' for which the deceleration parameter vanishes, as a 
function of the matter density in the DGP models, for the values of the matter density that 
are consistent (up to 99% confidence level) with type la supernova and CMB data (solid 
line), and the analogous quantity in the standard cosmological model (dashed line). In both 
cases fltot = 1- 



